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Cobalt-containing deposits from Cobalt carbonyl are experimentally produced and their composition is mea- 
sured. The Cobalt concentration is found to be readily tunable between 20 and 70 at.% by variation of the 
dwell time. The variations in metallic concentration are explained by co-deposition of hydrocarbons present 
in the chamber background pressure. 



I. INTRODUCTION 



II. EXPERIMENTAL METHODS 
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Focused electron-beam induced deposition (FEBID) 
is a promising technology for the maskless deposi- 
tion of metallic structures on a nano-scalei~— . The 
electron-beam is scanned in the presence of adsorbed 
organo-metallic precursors, leading to dissociation of the 
molecules and desorption of the volatile by-products, 
thus defining local, metal-containing structures in a range 
defined by the beam diameter and the primary electron 
interaction volume from which the emitted electrons, sec- 
ondaries and backscattered, arises. Deposition of vari- 
ous ferromagnetic metals has been investigated, such as 
Ironi~— , Nickel^ or Cobalt— . Clean magnetic metal de- 
posits were achieved using UHV setups 6 -, or high beam 
currents on relatively thick (300 nm) Si02 insulating 
films 1 —, or beam heating effects^^. 

However, it is in certain cases desirable to deposit a 
composite nanostructure having the magnetic metal em- 
bedded as nanocrystals in a carbonaceous matri x 10 i 13 
serving for mechanical stability and as oxidation barrier 
in ambient or liquid atmospheres, for instance in mag- 
netic scanning probe applications^. Furthermore, this 
composite material is of particular interest for the fab- 
rication of sub-micron Hall sensors, as this material ex- 
hibits a large extraordinary Hall effect (EHE)i^. A pre- 
cise control of the deposition process is thus required, 
as it dominantly influences the magnetic sensing prop- 
erties of the deposit^. Here, we present a study of 
the deposition process in planar deposits obtained from 
dicobalt-octacarbonyl (Co2(CO)s). While the deposition 
of cobalt-containing pillars, where the electron-beam is 
constantly illuminating a single spot on the substrate, 
has been investigated in depth ^ 11 ' 13 , only few studies 
of the linear or planar deposition of Cobalt, where the 
beam is raster scanned over a certain area, have been 
published so far~. We find metallic concentrations vary- 
ing between 20 and 70 at.% and explain these findings by 
the co-deposition of hydrocarbons present in the chamber 
background. 



a ) Electronic mail: laurent.bernau@empa.ch 



The FEBID experiments were carried out at room tem- 
perature in a Hitachi S-3600N scanning electron micro- 
scope (SEM), equipped with a tungsten filament elec- 
tron source. Silicon substrates covered by a 200nm-thick 
SiC>2 layer were used as a deposition substrate. Home- 
made modifications to the SEM allowed placing an in- 
ternal reservoir for the precursor, connected to a syringe 
nozzle pointed towards the substrate. Dicobalt octacar- 
bonyl (C02C08, CAS 10210-68-1) was used as a precur- 
sor molecule. In order to prevent spontaneous decom- 
position of cobalt carbonyl, the precursor is usually sta- 
bilized with 5-10 wt.% hexane (Sigma-Aldrich GmbH, 
Buchs, Switzerland). However, the hexane stabilization 
proves to be a rich source for carbon co-deposition dur- 
ing FEBID. When stored in an atmosphere of argon 
(VWR International, Dietikon, Switzerland), cobalt car- 
bonyl yields deposits of higher cobalt content when dis- 
sociated by FEBID. 

The precursor flux during the experiments was about 
4.4xl0 15 molecules/s, as estimated from mass loss mea- 
surements. According to our gas flow MC simulations 16 , 
this translates into 1.5 xlO 17 cm _2 s _1 impinging on the 
FEB irradiated spot. The residual chamber pressure was 
estimated from pressure measurements carried out with- 
out precursor, after degasing of the experimental setup. 
In the presence of an internal precursor reservoir, the 
monitored chamber pressure eventually saturates. De- 
position was carried out after stabilization of the cham- 
ber pressure, and the residual pressure was estimated to 
lxlO -5 mbar. 

We used the XENOS lithography system to control the 
electron beam during deposition. The deposition control 
parameters are the nominal dwell time i^, i.e. the time 
the beam irradiates a pixel before moving to the next, 
the inter-pixel distance, defined as the distance between 
to adjacent pixels irradiated by the beam, and the re- 
fresh time t r , which measures the time elapsed before 
the next iteration of the structure irradiation (see figure 
Dividing the total irradiation time by the total ex- 
posed surface yields the deposition dose in C/cm 2 . In our 
experiments, the total deposition dose was kept constant 
to 10C/cm 2 , while the dwell time was varied over two or- 
ders of magnitude (500ns to 50/is). The refresh time was 
always kept >10ms, to ensure full replenishment with 
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FIG. 1. Scanning parameters. The rectangular deposits are 
realized using a serpentine deposition path. Aa; and Ay de- 
note the increment in x and y direction and are identical 
throughout the experiments. Each pixel is irradiated for a 
dwell time td- The refresh time is the time elapsing between 
to subsequent irradiations of the same pixel and is given by 
the multiplication of the number of pixels with the dwell time. 

precursor molecules between two irradiation iterations. 
The inter-pixel distance was set to 30nm. This com- 
pares to a beam diameter (FWHM) estimated to 70nm 
(assessed experimentally using BeamMetrjii), for the ac- 
celeration voltage (25kV) and beam current (InA) used. 
This leads to an overlap of roughly 2, and the effective 
dwell time is the double of the nominal dwell time. The 
total irradiation time per pixel is then the double of the 
effective dwell time, as the overlap is present in both x- 
and y-direction. 

The deposits were characterized with regard to the 
composition using energy dispersive x-ray spectroscopy 
(EDX) , using a 3keV probe. The accuracy of the compo- 
sitional values obtained by this mean was checked to be 
correct within ±5at.% by calibration measurements on 
Cobalt-carbonate (C0CO3). The height of the deposits 
was assessed by atomic force microscopy (AFM). 



III. RESULTS 
A. Composition 

The metallic concentration of the deposits is found to 
be highly dependent on the dwell time for structures with 
a width of 600nm, as shown in figure [2l For dwell times 
below l/is, Cobalt concentrations as low as 20at.% are 
found, while dwell times of 20/is or more yield deposits 
with above 60at.% Cobalt. For deposits with a charac- 
teristic dimension of 3^m however, the Cobalt concentra- 
tion is within 60-70at.% independently of the dwell time 



FIG. 2. Composition of the deposits as measured by EDX for 
varying dwell times. Filled symbols are for (3/xm) 2 structures, 
open symbols are for (600nm) 2 deposits. 



used, with a tendency for lower metallic concentrations 
for higher dwell times. 

The balance of the composition contains Oxygen and 
Carbon. For the 600nm squares, the Carbon content 
scales inversely with the Cobalt content, decreasing from 
about 70at.% for dwell times of 500ns down to 22at.% for 
dwell times of 50/is, whereas the Oxygen level remains 
stable at around 20at.%. For the 3^tm squares, the Oxy- 
gen level remains stable within the measurement error 
at around 15at.% and a Carbon level of about 20at.% is 
found. 

The EDX measurements were performed by focusing 
a 3keV electron beam on the center of the deposit, so 
that the escape cone for the X-rays has a radius of well 
below 50nm. Thus, the discrepancy in the concentrations 
found for the 600nm squares cannot be attributed to edge 
effects during the measurement. 



B. Deposition rates 

The deposit heights are reported in figure |3] The cor- 
responding deposition rates (right-hand axis of figure [3]) 
are obtained by dividing the heights by the total irradi- 
ation time per pixel, as defined in the experimental sec- 
tion. As the total irradiation time per pixel, the beam 
current and the beam diameter are constant through- 
out all experimental series, the deposit heights and the 
corresponding deposition rates scale linearly. For the 
600x600nm 2 squares, the deposition rates are found to 
increase strongly for low dwell times (below 5/is), and are 
nearly constant for dwell times above this threshold. For 
the 3x3^tm 2 squares however, nearly constant deposition 
rates are found for all dwell times investigated. 

As a comparison, the contamination deposition rates, 
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FIG. 3. Heights of the investigated deposits, for (3/xm) 2 (full 
squares) and (600nm) 2 (full circles) Co2(CO)8 deposition. Er- 
ror bars show spread in height measurements along three 
adjacent AFM profiles. The contamination deposition was 
performed at 2xl0~ 5 mbar (upwards pointing triangles) and 
1 x 10 _5 mbar (downwards pointing triangles) . The right-hand 
axis gives the corresponding growth rates. 



TABLE I. Model parameters for Co2(CO)s and two hydro- 
carbons. Two molecules are taken as example for the fit 
of the contamination deposits: benzene (CsHg) is taken as 
an example of a "light" hydrocarbon, as its dissociation un- 
der electron exposure has been investigated in-^; 1,3-Bis(3- 
phenoxyphenoxy)benzene (C30H22O4) is a polyphenyl ether 
(PPE) used in the oil diffusion pump and is taken as an ex- 
ample of a "heavy" hydrocarbon. 
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a calculated from mass loss measurement and MC simulation of 

precursor distribution. 
b approximated from chamber background pressure. 
c estimated from dimensions of undissociated precursor molecule. 
d estimated from molar mass and density of undissociated 

molecule. 

estimated from molar mass and deposit density using the 
approach described inli. (The estimated decomposed volume of 
Co2Ca;0 B compares with 6.3xl0 -2 nm 3 used for the FEBID 



product of W(CO) 6 



obtained by depositing 600x600nm 2 squares in the ab- 
sence of any precursor, are shown, for background pres- 
sure of 1 and 2 x 10 _5 mbar. A strong dependency on the 
dwell time is found in this case, and the absolute value of 
the deposition rate at the same dwell time setting is found 
to depend linearly on the background chamber pressure, 
i.e. on the amount of hydrocarbon "precursor" available 
for the deposition process. 



IV. DISCUSSION 



The growth rate is described by the model^: 



R = 



t d 



n(t) dt 



(1) 



where V is the volume of the decomposed molecule, a is 
the decomposition cross-section, / is the PE flux and n(t) 
describes the surface coverage with precursor molecule at 
a given time. Disregarding surface diffusion, the precur- 
sor coverage is described by: 
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Assuming full replenishment between successive irradia- 
tions, the growth rate equation can be solved analytically 
by integration of equation [2] 



R = Vaf 



{n 



1 - exp (-k d t d ) 
k d t d 



(3) 



where n out and n„ represent the molecular coverage in- 
side and outside the electron beam, and k d is the time 



constant of the dissociative depletion process. Similarly, 
a time constant k r for the precursor replenishment be- 
tween irradiations can be introduced. 



k d = sJ/no + 1/t + af (4a) 

k r = s J /no + 1/t (4b) 

n in = sJ/ (sJ/n + 1/t + af) (4c) 

n ou t = sJ/ (sJ/n + 1/t) (4d) 

The condition for full replenishment is fulfilled for 
k r ■ t r 3> 1, i.e. when the refresh time is much higher 
than the time constant for the replenishment process. 
While the impinging flux of precursor molecules J can be 
estimated from mass loss measurements and MC simu- 
lations, the monolayer coverage uq approximated using 
the molecule's size, the sticking coefficient s and the res- 
idence time t are unknown molecular parameters in the 
case of Co2(CO)s- Using as a typical value T=lms and 
assuming s=l, the refresh time used in our experiments 
proves sufficient to fulfill full replenishment conditions. 
The growth rates for the 3/im 2 Cobalt carbonyl and the 
contamination deposits series were fitted using the values 
given in table U and the model curves are shown in figure 
S] and [5] (see table |TT| . Although the growth rate model 
used (equ. [3]) does not account explicitely for surface dif- 
fusion of adsorbed molecules, diffusive effect contribute 
by variations of the effective residence time. 

In figure [31 it can be seen that the growth rates for 
(600nm) 2 squares deposited from Co2(CO)g are discon- 
tinuous, presenting similarities with hydrocarbon disso- 
ciation for low dwell times (high dependence on dwell 
time), whereas the growth rates at higher dwell times 
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TABLE II. Rate fits for cross-section and residence time using 
the growth rate model. The fits are based on equ. [3] using 
the values given in tabled] 



precursor 

Co 2 (CO) 8 

PPE 

Benzene 



nm 

(.00495 ± .00078^ 
(2.3 ±3.0) 
(2.09 ± 2.26^, 



(LIS 

720 ± 150 
90 ±9 
160 ± 20 



19.7 
111 
62.5 



As a comparison, a measured a Benzene of 0.35nm 
:„18 



is reported 



b A fitted a Co2(CO)8 of 5x10 3 nm 2 was presented in 



(Utke,Purruker,MNE07). 
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FIG. 4. Growth rates of (3^m) 2 deposits from Co 2 (CO)s. 
The total dose was maintained constant at 10C/cm 2 . The 
solid line represents a fit of equation [3] The fit is extended 
to 60s and the growth rate from a steady-exposure pillar is 
indicated for comparison. 



approach the values in figure |U We propose that this be- 
haviour is due to co-deposition of both Cobalt carbonyl 
and background pressure hydrocarbons, which reach the 
dissociation spot by surface migration. For the 3/im 2 
deposits, surface diffusion can be neglected, as the grow- 
ing deposit represents a barrier for surface diffusion of 
hydrocarbon species. 

In the presence of two precursors, competitive adsorp- 
tion can be modelled using: 
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Solving this differential equation and assuming full re- 
plenishment (n out — limt^oo [n r (t)]) yields: 
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FIG. 5. Growth rates of (600nm) 2 contamination deposits at 
1 and 2xl0 _5 mbar. The solid line represents a fit of equation 
[3]to the values at 2xl0 _5 mbar using polyphenyl ether (PPE) 
or benzene as example hydrocarbon molecules. Using the 
fitted parameters, the model growth rates at lxl0 _5 mbar 
are drawn (dotted line). 



where n* represents the precursor surface coverage from 
the one-molecule model described by the constants in 
equation [4] Comparing with equation [4] the presence of 
the second precursor doesn't affect the dissociation and 
replenishment time constants in the model. However, 
the surface site occupancies are limited by the presence 
of the competing precursor species. The growth rate is 
then expressed as the sum of the growth rate contribution 
for each species: 
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id 
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n 1 {t)dt + V 2 a 2 ft- 1 / n 2 {t)dt{7) 



The integrals were solved numerically and the result 
is shown in figure [SI The model predicts reasonably 
the change in growth rate observed for the (600nm) 2 se- 
ries. Figure [7] shows the variation of dissociation yield 

id 

(Yi — <Jit J rii(t)dt) per electron for the Cobalt car- 
o 

bonyl and hydrocarbon species. As the dwell time in- 
creases, the adsorbed hydrocarbons are locally depleted, 
as they are dissociated efficiently, and this increases the 
relative contribution of Cobalt carbonyl to the deposit 
growth. This is similar to the switching between elec- 
tron beam induced etching and deposition described in2£, 
where selective depletion was achieved by varying the 
beam current used. 

Under electron irradiation, the carbonyl groups of the 
Co2(CO)s are partially dissociated from the metallic 
atoms and escape as volatile molecules, or are themselves 
dissociated under electron irradiation and get fixed in 
the growing deposit 1 . Thus, prediction of the deposit's 
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FIG. 6. Deposition rates for (600nm) 2 squares. The growth 
rates for (3/im) 2 structures and for contamination structures 
are given for comparison. The line is the model prediction as 
described by equation [7] 



FIG. 8. Composition of 600nm squares and model prediction. 
Full line: PPE, dotted line: benzene. 
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FIG. 7. Calculated yields for Cobalt carbonyl, benzene and 
PPE as example hydrocarbons. 



composition is difficult, as the dissociation pathways are 
manifold, resulting in a complex process. However, the 
dissociation product of Cobalt carbonyl in our experi- 
mental conditions can be estimated using the composi- 
tion of the (3/zm) 2 deposits shown in figure [21 Based on 
these EDX measurements, we approximate the dissocia- 
tion product of Cobalt carbonyl in our experimental con- 
ditions to Co2C0.6iO0.46- Furthermore, we assume that 
carbon contamination is dissociatively deposited without 
volatile carbon-containing fragments escaping. It is then 
possible to model the composition of the deposit where 
co-deposition of hydrocarbon species occurs, using: 



where Y denotes the dissociation yield, the lower term 
represents the number of deposited atoms per electron, 
and the upper term the number of deposited Co atoms 
per electron. The result is shown in figure [8] and com- 
pared with the composition found for the (600nm) 2 de- 
posits. The two-precursor model describes the variations 
in composition quite well. As depletion of hydrocarbons 
is more pronounced, longer dwell times lead to cobalt-rich 
deposits. For even longer dwell times, a slight decrease of 
Cobalt content is predicted, as depletion of Cobalt car- 
bonyl occurs. 

I n?! 1 ?! 12 ! 1 ? , Co concentrations between 35at.% and 
90at.% are reported. The differences are usually at- 
tributed to thermal effects by local electron-beam in- 
duced heating. However, for low beam currents (< InA) 
and Si02 as a substrate, MC simulation of electron-beam 
heating effects predict negligible local heating. On the 
other hand, the proposed co-deposition model explains a 
large scatter in metallic composition even at background 
pressures of 10~ 6 mbar, lower than the one investigated 
here, depending on the ratio of adsorbed hydrocarbon 
molecules participating in the deposition. 



V. CONCLUSION 

The importance of contamination incorporation was 
shown for tungsten from WF6 in 21 -, where plasma- 
cleaning of the surface and removal of surface-adsorbed 
hydrocarbons would lead to deposits of higher purity, as 
well as ir£ for the organic precursor Fe(CO)s, where an 
ultra-high- vacuum environment increased the deposit pu- 
rity. The question is discussed in the topical review in^ 2 -. 
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Here, it was shown how to take advantage of the con- 
comitant deposition of residual hydrocarbons in order to 
tune the Cobalt content of planar deposits from Cobalt 
carbonyl within a range of 20 to 70 at.%. The findings 
are especially relevant for using such deposits as magnetic 
sensing material, as the Hall resistivity depends strongly 
on the metallic content of Cobalt deposits, with a strong 
peak at around 65at.%22. 
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